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Abstract 



^\ . In this paper, we consider the short time strong solution to a simphfied hydrodynamic flow 

04 ' modeling the compressible, nematic liquid crystal materials in dimension three. We stablish a 

criterion for possible breakdown of such solutions at finite time in terms of the temporal integral 
p ^ ' of both the maximum norm of the deformation tensor of velocity gradient and the square of 

. maximum norm of gradient of liquid crystal director field. 

^' _ 

■ 1 Introduction 

Nematic liquid crystals are aggregates of molecules which possess same orientational order and are 
made of elongated, rod-like molecules. The continuum theory of liquid crystals was developed by 
^ I Ericksen [9] and Leslie [28] during the period of 1958 through 1968, see also the book by de Gennes 

' [TT]. Since then there have been remarkable research developments in liquid crystals from both 

theoretical and applied aspects. When the fluid containing nematic liquid crystal materials is at 
ly^ ', rest, we have the well-known Ossen- Frank theory for static nematic liquid crystals. The readers 

I can refer to the poineering work by Hardt-Lin-Kinderlehrer [12] on the analysis of energy minimal 

CD ' configurations of namatic liquid crystals. In general, the motion of fluid always takes place. The 

so-called Ericksen-Leslie system is a macroscopic continuum description of the time evolution of 
the materials under the influence of both the flow velocity field u and the macroscopic description 
of the microscopic orientation configurations d of rod-like liquid crystals. 
^ . When the fluid is an incompressible, viscous fluid, Lin [18] first derived a simplified Ericksen- 

I Leslie equation modeling liquid crystal flows in 1989. Subsequently, Lin and Liu [191 ED] made 

some important analytic studies, such as the existence of weak and strong solutions and the partial 
regularity of suitable solutions, of the simplified Ericksen-Leslie system, under the assumption that 
the liquid crystal director field is of varying length by Leslie's terminology or variable degree of 
orientation by Ericksen's terminology. 

When the fluid is allowed to be compressible, the Ericksen-Leslie system becomes more com- 
plicate and there seems very few analytic works available yet. We would like to mention that very 
recently, there have been both modeling study, see Morro [29j, and numerical study, see Zakharov- 
Vakulenko [36j . on the hydrodynamics of compressible nematic liquid crystals under the influence 
of temperature gradient or electromagnetic forces. 

The main aim of this paper, and the companion paper [T7] as well, is an attemption to initiate 
some analytic study for the flow of compressible nematic liquid crystals. We will mainly address 
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several issues on the strong solutions. More precisely, we will focus on the blow-up criterion on 
strong solutions in this paper. 

Now we start to describe the problem. Let Q C M'^ be either a bounded smooth domain or 
the entire M^, we will consider a simplified version of Ericksen-Leslie equation that models the 
hydrodynamic flow of compressible, nematic liquid crystals in fi: 

pt + V-(pn) = 0, (1.1) 
put +pu-Vu + V(P(p)) = Cu-Ad- Vd, (1.2) 
dt + u-Vd = Ad+\Vd\'^d, (1.3) 

where p : Q, ^ M+ is the density of the fluid, ti : — )• is the fluid velocity field, P{p) : — s- 
denotes the pressure of the fluid, d : — )• S*^ represents the macroscopic average of the nematic 
liquid crystal orientation field, V • (= div) denotes the divergence operator on R^, and C denotes 
the Lame operator defined by 

Cu = pAu + {p + \)\/{V -u), (1.4) 

where p and A are the shear viscosity and the bulk viscosity coefficients of the fluid repsectively, 
which are assumed to satisfy the following physical condition: 

p>0, 3A + 2/i>0. (1.5) 

The pressure P{p), as a given continuous function of p, is usually determined by the equation of 
states. Through this paper, we assume that 

P : [0, +00) — )• M is a locally Lipschitz continuous function. (1-6) 

Notice that (jl.ip is the equation for conservation of mass, (jl.2p is the linear momemtum equation, 
and (11. Sp is the angular momentum equation. We would like to point out that the system (jl.ip - (jl.3p 
includes several important equations as special cases: 

(i) When p is constant, the equation (jl.ip reduces to the incompressibility condition of the fluid 
(V • u = 0), and the system (ll.ip - ljl.Sp becomes the equation of incompressible flow of namatic 
liquid crystals provided that P is a unknown pressure function. This was previously proposed by 
Lin |18j as a simplified Ericksen-Leslie equation modeling incompressible liquid crystal fiows. 

(ii) When d is a constant vector field, the system ()l.ip - ()1.2p becomes a compressible Navier- 
Stokes equation, which is an extremely important equation to describe motion of compressible 
fluids. It has attracted great interests among many analysts and there have been many important 
developments (see, for example. Lions [26], Feireisl |10j and references therein). 

(iii) When both p and d are constants, the system (jl.ip - (jl.2p becomes the incompressible Naiver- 
Stokes equation provided that P is a unknown pressure function, the fundamental equation to 
describe Newtonian fluids (see. Lions [25J and Temam [30 1 for survey of important developments). 

(iv) When p is constant and u = 0, the system (jl.ip - (|1.3p reduces to the equation for heat flow 
of harmonic maps into S"^ . There have been extensive studies on the heat flow of harmonic maps in 
the past few decades (see, for example, the monograph by Lin- Wang [23] and references therein). 

Prom the viewpoint of partial differential equations, the system ()l.ip - ()1.3p is a highly nonlinear 
system coupling between hyperbolic equations and parabolic equations. It is very challenging to 
understand and analyze such a system, especially when the density function p may vanish or the 
fluid takes vacuum states. 
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In this paper, the system (|l.ip - (ll.3p wih be studied along with the initial condition: 

{p,u,d) = {po,uo,do), (1.7) 

and one of the following three types of boundary conditions: 

(i) Cauchy problem: 

= M^, p,u vanish and d is constant at infinity (in some weak sense). (1-8) 

(ii) Dirichlet and Neumann boundary condition for {u,d): 17 C is a bounded smooth domain, 
and 

where i' is the unit outward normal vector of d^l. 

(iii) Navier-slip and Neumann boundary condition for {u,d): (7 C is a bounded, simply con- 
nected, smooth domain, and 



(u ■ v, curlii X u, -—] 
ou 



where curlu = V x u denotes the vorticity field of the fluid. 

In order to state the definition of strong solutions to the initial and boundary value problem 
prT]) - (fL3|) . together with ([THD or or pTTO]) . we introduce some notations. 

We denote 

f dx = / / dx. 



For 1 < r < oo, denote the spaces and the standard Sobolev spaces as follows: 

L'^ = L^{n), D^'" = {-u G Ll^{n) : ||V^u||l'- < oo} , 

DI = ^ue : ||Vn||i2 < oo, and satisfies (fL8]) or (fL9l) or (ll.lOj) for the part of n|, 

Hl = L^r\Dl, \\u\\^k,r = \\V^u\\Lr. 



Denote 
and let 



QT = nx [0,r] (T>0), 



V{u) = \{Vu + {Vuf) 



denote the deformation tensor, which is the symmetric part of the velocity gradient. 

Definition 1.1 For T > 0, {p,u,d) is called a strong solution to the compressible nematic liquid 
crystal flow ^.l\) - fT73\} in U x (0,T], if for some q G (3,6], 

< p e C([0, T]; H^i'^ n H^), Pt G ^([0, T]; Q L^); 
«GC([0,r];Z)2nDi)n^'(0,r;Z)2,g)^ ntGL2(0,T;Di), € L-(0, T; L^); 

Vd G C{% T];H^) n L\0, T; H^), dt G C{% T]-H^)\^ L^{0, T; H^), \d\ = 1 in Q^, 
and {p,u,d) satisfies lil. 1\} - [T7^} a.e. m x (0,T]. 
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For the existence of local strong solutions associated with the three types of boundary conditions, 
we have obtained the following theorem in the paper |17j . 

Theorem 1.2 Assume that P satisfies lll.6\) . po > 0, po G W^''^ f] f] for some q G (3,6], 
uq G D^PIDq, Vdo G and \do\ = 1 inQ. If, in additions, the following compatibility condition 

Cuq - V{P{po)) - Ado ■ Vdo = Vpog for some g G M^) (1.11) 

holds, then there exist a positive time Tq > and a unique strong solution {p,u,d) of ^.l\) - fT73\) . 
[LTD together with fTgj or or [rW\) in n x {0,To]. 

We would like to point out that an analogous existence theorem of local strong solutions to 
the isentropic compressible Naiver-Stokes equation, under the first two boundary conditions (jl.Sp 
and (jl.9p . has been previously established by Choe-Kim [3] and Cho-Choe-Kim [4J. A byproduct 
of our theorem 11.21 also yields the existence of local strong solutions to the isentropic compressible 
Navier-Stokes equation under the Navier-slip boundary condition (jl.lOp . 

In dimension one, Ding-Lin- Wang- Wen [7j have proven that the local strong solution to (jl.ip - 
(jl.Sp under (|1.7p and (|1.9p is global. For dimensions at least two, it is reasonable to believe that 
the local strong solution to (jl.ip - (jl.3p may cease to exist globally. In fact, there exist finite time 
singularities of the (transported) heat flow of harmonic maps (jl.Sp in dimensions two or higher (we 
refer the interested readers to [23j for the exact references). An important question to ask would 
be what is the main mechanism of possible break down of local strong (or smooth) solutions. 

Such a question has been studied for the incompressible Euler equation or the Navier-Stokes 
equation by Beale-Kato-Majda in their poineering work [1], which showed that the L°°-bound of 
vorticity V x u must blow up. Later, Ponce [27] rephrased the BKM-criterion in terms of the 
deformation tensor P(u). 

When dealing with the isentropi^ compressible Navier-Stokes equation, there have recently 
been several very interesting works on the blow up criterion. For example, if < T,,, < +oo is the 
maximum time for strong solution, then (i) Huang-Li-Xin [H] established a Serrin type criterion: 
lim^^r^ (l|divu||ii(o,T;L°°) + II/'^'"IIl»(o,T;L'-)) = oo for I + ^ < 1, 3 < r < oo; (ii) Sun- Wang-Zhang 
[HI] , and independently [H], showed that if 7p > A, then lim7--|-j'^ IIpIIl°°{o,T;L°°) = co; and (iii) 
Huang-Li-Xin [15] showed that limx^x* II^(^)IIli(o,T;L°°) = oo. 

When dealing the heat flow of harmonic maps (|1.3p (with n = 0), Wang [32] obtained a Serrin 
type regularity theorem, which implies that if < < +oo is the first singular time for local 
smooth solutions, then lirxiTfTt ||V(i||2,2(o,T;L°°) = c». 

When dealing with the incompressible nematic liquid crystal flow, Lin-Lin- Wang [24j and Lin- 
Wang [22] have established the global existence of a unique "almost strong" solutioro for the initial- 
boundary value problem in bounded domains in dimension two, see also Hong [13j and Xu- Zhang 
|34j for some related works. In dimension three, for the incompressible nematic liquid crystal flow 
Huang- Wang |16j have obtained a BKM type blow-up criterion very recently, while the existence 
of global weak solutions still remains to be a largely open question. 

Motivated by these works on the blow up criterion of local strong solutions to the Navier-Stokes 
equation and the incompressible nematic liquid crystal flow, we will establish in this paper the 
following blow-up criterion of breakdown of local strong solutions in finite time. 



^namely, P{p) — ap"' for some a > and 7 > 1. 

^that has at most finitely many possible singular time. 
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Theorem 1.3 Let (p,u,d) be a strong solution of the initial boundary problem ji. ip -i TO]) . 7| j 
together with /il.8\) or \1.9\) or \1.10\) . Assume that P satisfies \1.6\) . and the initial data {pQ,UQ,do) 
satisfies //O < < +00 is the maximum time of existence, then 

r (llP(n)IUco + llVdIlioo) dt = 00. (1.12) 
^0 

We would like to make a few comments now. 

Remark 1.4 (a) In [17], we also obtained a blow-up criterion of (|l.ip - (jl.3p under the initial con- 
dition (II. 7p and (II. Sh or (II. 9p in terms of p and Vd. Namely, if 7p > 9A and < < +00 is the 
maximum time of existence of strong solutions, then 

1™ (|I/5|Il°°(o,T;L-) + l|V(i||z,3(o,r;L°°)) = +00. 

(b) For compressible liquid crystal flows without the nematicity constraint {\d\ = ijl, Liu-Liu |21j 
have recently obtained a Serrin type criterion on the blow-up of strong solutions. 

(c) It is a very interesting question to ask whether there exists a global weak solution to the 
initial-boundary value problem of ()l.ip - ()1.3p in dimensions at least two. In dimension one, such an 
existence has been obtained by Ding- Wang- Wen [8]. 

We conclude this section by introducing the main ideas of the proof, some of which are inspired 
by some of the arguments on the isentropic compressible Navier-Stokes equation by |15j and [31j . 

(1) It is well-known that the bound of 11^(^)1 1^,1 2,00 yields that ||/o||Lj°°Lg° is bounded from the equa- 
tion (11. ip . See Lemma 2.1. 

(2) We observe that in the equation (jl.Sp . the bound of (||^(^i)||LiLcx) + ||V(i||^2^oo) yields that 
II V(i||Lj«Lr; is bounded for any 2 < r < +00, which is a crucial ingredient to obtain higher order 
estimates of p, u, d. See Lemma 2.3. 

(3) Due to the possible vacuum state of p, the strong nonlinearities of the convection term u ■ Vu 
and the induced stress tensor Ad • Vd, in order to obtain control of (||Vp||ioo^2 + ||Vu||^oo^2 + 
II V^(i||j;^cx)^2 ), we estimate (||\/p'w||l2j;^2 + ||V(it ||j;^22^2 ) by combining an energy estimate of the equa- 

^^^^^^ ^ £ 3^ i cc 

tion (jl.2p in terms of the material derivative u = ut + u- \/u with second order energy estimates of 
both ([O]) and ([O]). See Lemma 2.4. 

(4) We estimate (|| V^u||2,^l2 + || V^d||2,^2,2 ) by combining thrid order estimate estimates of (jl.2p 
and (jl.3p with ff^-estimate of the Lame equation and i7^-estimate of the harmonic map equation. 
See Lemma 2.6. 

(5) Finally, we obtain the estimate of ||V/?||^oo^9 for 3 < g < 6 in terms of ||u||^2£)2,g. To do 
it, we employ the elliptic estimate of the equation satisfied by the effective viscous flux G = 
{2p + A)divti — P{p) and the bound of ||V^(i||2,2^2 and || Vnt||j;^22^2 . See Lemma 2.7. 

We would like to point out that during all these arguments, specific forms of the pressume 
function P play no roles, it is the local Lipschitz regularity of P that is relevant. 

Acknowledgement. The first two authors are partially supported by NSF grant 1000115. The 
work was completed during the third author's visit to the University of Kentucky, which is partially 
supported by the second author's NSF grant 0601162. The third author would like to thank the 
department of Mathematics for its hospitality. 

^the right hand side of equation (11.31) is replaced by Ad + f{d) for some smooth function f : R'' — >■ R'^, e.g. 
/(d) = (Idp - l)d. 
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2 Proof of Theorem D 



Let < < oo be the maximum time for the existence of strong solution (p,u,d) to (ll.ip -( fL3]) . 
Namely, {p,u,d) is a strong solution to (jl.ip - (jl.3p in x (0, T] for any < T < T*, but not a 
strong solution in 17 x (0,T*]. Suppose that (I1.12P were false, i.e. 



Mo := / ' (||:P(n)||L- + llVdIlioo) dt < oo. (2.1) 
Jo 



The goal is to show that under the assumption (12. ip . there is a bound C > depending only on 
Mq, po, UQ,do, and T* such that 



sup 

0<i<T* 



max(||/9||p^i,r + ||/5t||Lr) + {\\y/put\\L2 + ||Vn||j^i) + (||dt||//i + ||Vd||jf2) 

r=2,g 



< C, (2.2) 



and 



/ * ihtWli + II«IId2,, + ||dt||^2 + llVdll^a) < C. (2.3) 
Jo 

With ()2.2p and p.3p . we can then show without much difficulty that is not the maximum time, 
which is the desired contradiction. 

Throughout the rest of the paper, we denote by C a generic constant depending only on po, uq, 
do, T^, Mo, A, p, and P. We denote by 

if there exists a generic constant C such that A < CB. For two 3x3 matrices M = (Mij), N = (Nij), 
denote the scalar product between M and by 

3 

For d : — 7- S^, denote by Vd ® Vd as the 3x3 matrix given by 

{Vd(g)Vd)ij = {Vid,Vjd), 1 < i,i < 3. 
The proof is divided into several steps, and we proceed as follows. 

Step 1. We will first establish L°°-control of p. More precisely, we have 

Lemma 2.1 Let < T* < +oo be the maximum time for a strong solution {p,u,d) to lil.l\) - fT7^) . 
[LID together with [TB^ or or iLm . If KJl]) and hold, then 

sup IIpIIloo < C. (2.4) 

0<t<T, 

Proof. This estimate is a well-known fact of (jl.ip and was proved by Huang-Li- Xin [15J (Lemma 
2.1). For the convenience of reader, we sketch it here. For any 1 < r < +oo, multiplying (II. ip by 
rp"^^^ and integrating over fi, we obtain 



d_ 
It 



J p^ dx = — J {u ■ V{p^) + rp'^divii) dx 

= — j (div (up^) + (r — l)/9'"divti) dx < (r — l)||divti||j;,oo j p' dx. 
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Thus 

d r — 1 

37 P L'- < divn LOO p 

at r 

This, (12. ip . Lemma 2.1, together with Gronwall's inequahty, imply 

sup < II/OoIIl'- exp ( / ||divu||ioodt ) < C, 

which, after sending r to oo, implies ()2.4p . This completes the proof. □ 

Step 2. We next establish the global energy inequality for strong solutions, namely. 

Lemma 2.2 Let < < +oo be the maximum time for a strong solution {p,u,d) to ^.l\) - fT73\) . 
( fTTP together with or or ^A^). If [JJl]) and hold, then for any <t <T^, the 

following inequality holds: 



< c 

Furthermore, we have 



{p\uf + \Vdf) {t)dx + j (\Vu\'^ + \M+\Vd\'^d\^^dxds 
/ m^^ol^ + iVdoP) dx + l . 



(2.5) 



V'^d\'^dxdt<C. (2.6) 

Jn 



Proof. Without loss of generality, we may assume that P(0) = 0. Since P is locally Lipschitz by 
()1.6p . it follows that P' is locally bounded on [0, +oo). Since p is bounded in $7 x [0,T,,) by ()2.4p . 
we then have that, on $7 x [0,T^,), 

\P{p)\ < \\P'{p)\\l^P<Cp{<C) (2.7) 
|V(P(p))| < ||P'(p)||^^|Vp| <C|V/,|. 

Since (jl.Sp and (jl.9p are easier to handl^, we outline the proof for the boundary condition (jl.lOp . 
Multiplying (|1.2p by u and integrating over ^2, we have 

^ j p{dt\u\'^ + u ■ V|np) dx + j {p\V x up + {2p + A)|divup) 

= j P{p)divudx - j u-Vd- Addx. (2.8) 

Here we have used the fact that Au = Vdivn — V x curl n, and the Navier-slip boundary condition 
pTTO]) to obtain 

J Cu ■ udx = j [{2p + A) V(divn) • u — pV x curl u ■ u]dx 
= — J {2p + X)\diYu\'^ + plcmlul"^) dx. 



^in fact, with respect to the boundary condition p.8p and (|1.9p . by integration by parts one has 
{^lAu + (a* + A)V(divu)) - udx = / (m|Vu|^ + (p + A)|divM|^) dx > \\/uf dx. 
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Hereafter we repeatedly use the following identity: 

y (V X n, curl u) dx = j {u, V x (curl u)) dx, Vn with curl u x = on dQ. 
By the formula of transportion, we have 

dt J 

By ()2.7p and Cauchy's inequality, we have 

I J P{p)divudx\ < j \P{p)\\divu\dx < J p\Vu\dx 

< e j |Vn|2 dx + (7(e) j < e j |Vn|2 dx + C(e), 
where we have used (j2.4p and the conservation of mass to assure 

j - \\p\\lA\p\\l°° < C. 
Putting these inequalities into (|2.8p . we obtain 

J p\u\'^dx + j\p\Vxu\^ + {2p + X)\diYu\^)dx<- j u-Vd-Addx + e J |Vnp (ix + C(e). (2.9) 



d_ 
di 



Since is assumed to be simply connected for the boundary condition (ll.lOp . we have the following 
estimate (see [33] for its proof): 

l|Vii||L2 ^ l|V X u\\l2 + ||divn||i2 Vn G H^{n) with u • i/ = on 90. (2.10) 

This, combined with (jl.Sp . implies that 

j {p\V X up + {2p + X)\divu\'^)dx > I y"(|V X up + |divup) dx > ^ J |Vup dx. (2.11) 

Thus, by choosing e = ^I'^-^P implies 

^ y /9|up dx + ^ J |Vtip dx<-J u-Vd- Addx + C. (2.12) 

Now, multiplying (II. 3p by (Ad + | Vdpd), integrating over Q, and using |^ = on dQ, we have 

~ J |Vdpdx + J |Ad+ |Vdpd|^dx = J u-Vd-Addx, (2.13) 

where we have used the fact that |d| = 1 in and hence 

j\dt + u- Vd, |Vdpd) dx = 0. 

Adding ([232]) and (l2T3]l together yields (f23]l . 



To see <^M), observe that ([21]) implies J^* \\Vd\\l^ dt < Mq so that 
/ * / \Vd\^dxdt < Mo - ( sup / \Vd\^dx] 

Jo Jn \ 0<i<T, J I 



< CM, 



1 + / iPoWol^ + iVdoP) dx 





where we have used ()2.5p in the last step. This and (I2.5P then imply 

r-T, r rT, r „ rT* 



Jn 



\Ad\^dxdt = [ ' [ \Ad+\Vd\^df dxdt+ [ * [ \Vd\ 
Jo Jn Jo Jn 



^ dx dt 



< CMo 



1+ / (/Ool^toP + iVdoH dx 



Since ^ = on dQ, the standard L^-estimate yields 

J iV^dp dx<C j (I Adl^ + I Vdp) dx. 
Thus (j2.6p follows easily, and the proof is complete. □ 

Step 3. We will establish L^L^-control of Vd for any 2 < r < +oo, a key ingredient for the higher 
order estimates of u, Vd. More precisely, we have 

Lemma 2.3 Let < < +oo be the maximum time for a strong solution {p,u,d) to il.l\) - [T73\) . 
(flTTj ) together with (fTgj or / fO) or / flTf^) . // / fTTTF]) ond (EJP /loW, i/ien /or any 2 < r < +oo, 
there exists a C > depending on MQ,uo,do,Q,n, and r such that 

sup llVdll^,. + / * / \Vd\''-^\V^d\'^dxdt<C. (2.14) 
o<t<T, Jo Jn 

Proof. Here we only consider the Navier-slip boundary condition (|1.10p . since the argument to 
deal the first two boundary conditions (jl.Sp and (|1.9p is similar and easier. Differentiating the 
equation (11. 3p with rrespect to x, we have 



Vdt - VAd = -V(n • Vd) + V(|Vd|^d). (2.15) 
Multiplying (I2.15P by r|Vd|''~^Vd and integrating over Q, we obtain 

|Vdrdx + r j (|Vdr-2|v2^|2^(^_2)|vd|'-2|v|Vd||2) dx 

--r j V(|Vd|2d)|Vdr"2vddx-r j Viu ■Vd)\Vd\'^'^Vddx 16) 



d_ 
'dt 



r 

+ 2 



3 

/ |Vdr-2(V(|Vd|2),i.)da = V/^ 
Jdn 



We can estimate (i = 1,2,3) separately as follows. For /i, since 

V(| Vdpd) = iVdpVd + V(| Vdp)d and d • Vd = 0, 
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we have 

For I2, we have 
h = - 



I^=r J \Vd\''+^dx <\\Vd\\l^ J iV^rdx. 



\Vd\''-'^ViU^Vjd, Vid) dx- j u- V(|V(ir) dx 

= -rj \Vd\''-^V{u) ■.Vd0Vddx + j {divu)\Vd\'' dx < \\V{u)\\l-- J \Vd\'' dx. 

The estimate of the boundary integral I3 is more dehcate. Let Iqq denote the second fundamental 
form of dQ: for any x G dQ, 

lanix){U, V) = -Vv{x){U, V), V[7, V e T,{dn). 

Let Vt denote the tangential derivative on dO.. Since |^ = (Vd, z^) = on Of], we have VT'(f^) = 

on dVt. Hence we have, on dVt, 

1 c)di 

-(V(|Vd|2),z^) = Vd-V{Vd,v)-Vv{Vd,Vd) = VTd-VT{-Q^)-Vi^C^Td,VTd) = lanlVrd, Vtc^). 



Therefore we have 



h=r [ \Vdr%n{yTd,VTd)da< [ {Vdl"- da. 
JdQ JdQ 



Ian Jan 
Applying the trace formula W^'^{Q,) C L^{dQ,) and Holder's inequality, we obtain 

^3<iiivdriiH/M< I \vd\^dx+ 1 \vdr'\v^d\dx 
<c J |vdr'dx+ J J \vdr^\v^d\^dx. 

Putting all these estimates into (j2.16p . we obtain 

^ j \Vd\''dx + ^ j \Vdr''\V^d\^dx<{\\Vd\\l^ + \\Viu)\\Loo+l) j \Vd\'dx. 
By Gronwall's inequality and (j2.ip . we obtain that for any < t < T^,, 

/ |Vd(t)rdx+ r [ \Vd\''-^\V'^d\'^ dxds 
J Jo Jn 

< j iVdol'-dx -6x^(^1^ \l + \\Vd\\lo. + \\V{u)\\L^)d?j < C. 
This completes the proof. □ 

Step 4. Estimates of {Vu,Vp,V^d) in LfLl{Q. x [0,T*]). First, for any function / on x (0,r*), 
let 

f = ft + u-Vf 
denote the material derivative of /. Then we have 
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Lemma 2.4 Let < < +00 he the maximum time for a strong solution {p,u,d) to il.l\) - nT3\) . 
[n\ j together with (Eg) or or WJUO . If n~Jl\) and (KM> hold, then 

sup (llV-ull^^ + ||V/)|||2 + llV^dll^^) + [ * [ + \Vdt\'^) dxdt < C. (2.17) 

o<t<T* Jo Jn 



Proof. To make the presentation shorter, here we only consider the difficult case: the Navier-slip 
boundary condition (jl.lOp . To obtain the estimates of u, we adapt some arguments by [15] (Lemma 
2.2). Multiplying ()1.2p by ii and integrating over fl, we obtain, 



/9|n|^ dx — J {Cu, ut) dx 

= Vn, /:u)dx- j u-Vu- V{P{p)) dx - j Uf V{P{p)) dx (2.18) 

- j u-Vu - {Ad, Vd) dx - j Uf {Ad, Vd) dx. 

Similar to the proof of Lemma 2.1, we have 

— J {Cu, Ut) dx = J [/i(V X curlu, ut) — (2^ + A)(V(divu), n^)] dx 

[/i(V X u, V X Uf) + {2fj, + A)(divn)(divnt)] dx 

j X up + (2/i + A)(divn)^] dx, 

where we have used the fact that ut ■ v = curln x = on dVt during the integration by parts. 
The terms on the right hand side of ()2.18p can be estimated as follows. 

j {u-Vu,Cu)dx = - n J{u-Vu,V X curl u) + (2/i + A) J {u ■ Vu,V{divu)) dx. (2.19) 
For the first term in the right hand side of ()2.19p , by using curl u x ly = on dil. and the formula 

u X curln = -V(|np) — u ■ Vu, 



1 d 
'2di 



we have 



— p y ■ ^ ^ ^^^'^ u) dx = —p J (curl u, V X (u ■ Vu)) dx = fi J curl n • V x (u x curl u) dx 

=fi J (curl u, (curl u ■ V)ti — (n • V)curl u + div (curl u)u — (div ti)curl u) dx 

=fi J ^P(n) : curln (8) curln — -div ti(curl?x)^^ dx < ||'D(n)||Loo ||Vn||^2, 

where we have also used the formulas 

V X (a X 6) = (6 • V)a — (a • V)6 + (div6)a — (diva)6, and div(curln) = 0. 
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To estimate the second term in the right hand side of (j2.19p . denote by n'^ = n — (li • i/)z/ the 
tangential component of u on dVt. Note that (ll.lOp imphes u = on dQ. Hence we have 

/(« • Vu, V(divn)) dx = / {{u ■ V)u,u)divuda - / ((Vn) : (Vti)*divu + • V((divu)2)) dx 

= / u"^ -Vriu ■ iy)diYuda — / Vi'{u'^ ,u^)diY uda 
Jdn Jdn 

1 -i 
[(Vn) : (Vn)Mivn - -(divu)'^] dx 

Iduiu^ ,u^)divuda - / [(Vn) : (Vu)*divu- ^(divu,)^] 
^Il^lli4(af7)l|divn||i2(gf^) + \\V{u)\\loc\\Vu\\12. 
By the trace formula H^{^) C U'{d^) for r = 2,4, the Poincare inequality (see [35]): 

IKIIl^ ^ II^'^^IIl^, Vu G i?"*" with li • 1/ = on 9$^, 
and Holder's inequality, we have 

ll«lli4(c»n)l|divn||L2(aQ) <||u||^i ||Vu||hi < ||Vu||^2(||Vn||L2 + ||V\||l2) 

<C(e)(l + ||Vn||i2) + e||V\||i2 
for small e > to be determined later. Thus we obtain 

j {u ■ Vu,Cu) dx < €\\V^u\\l2 + C(e)(l + ||Vn||^2) + C\\V{u)\\loo \\Vu\\l2. (2.20) 
The remaining terms in the right hand side of ()2.18p can be estimated as follows. 
- j u-Vu-V{P{p))dx 

P{p){{u ■ V)n, iy)da+ / {P{p)u ■ V(divn) + P(p)(Vn) : (Vuf ) dx 



an 

t 



[P{p){u • V)(n • I/) - P{p)l9n{u^, u^)] da + / P{p){Vu) : {Vuf dx 
on J 

+ / P{p){u-iy)dwuda- [V{P{p)) ■ udiv u + P{p){dw uf)] dx 
Jdn J 

P{p)l9n{u\u^) da + [[P{p)i{Vu) : (Vn)* - (divn)^)] dx (2.21) 



V(P(p)) • ndivndx 

^\Mh{dn) + l|Vn||i2 + ^ |VH|n||divn| dx 

<||Vn||i2 + ||P(n)||i3||n||i6||Vp||,.2 

<||Vn||i2 + ||P(n)|||^||Vn|||,||VHU2 

<1 + (llP(n)IU^ + 1) ||Vn||i2 + ||Vn||i2||Vp||i2, 

where we have used (|2.7p and the Sobolev and Poincare inequalities (see p5l): 

\\u\\L<i ^ i\\u\\L2 + l|Vn||i2) < ||Vn||i2, Vn G H'^^Q) with u- = on dQ. 
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Since ()2.4p and (jl.ip also imply 

mp))t\ < \\P'{p)\\L^\pt\ < \\P'ip)\\L^ip\Vu\ + \Vp\\u\) < Ci\Vu\ + \u\\Vp\), (2.22) 
we have 

— J ut ■ VP{p) dx = — J P{p)div udx — j {P{p))td\\ udx 
-Jt j P(p)^''''^dx + C ji\u\\Vp\\dwu\ + \Vu\'^)dx 

<^ J Pip)diyudx + C[\\Vu\\l, + \\Viu)h3\\uh4'^pU2] (2.23) 



d f 

<— P{p)divudx + C 



\\Vu\\i, + \\V{u)\\l^\\Vu\\l,\\Vph2 



1 /" 

<- J P(p)divndx + C[l + (l + ||P(^x)||L^)||Vn||i. + ||V^x||i.||Vp||i.], 
where we have used the Poincare inequality and Holder's inequality. 

-J u-Vu- {Ad,Vd)dx <\\u\\L6\\Vu\\L(i\\Ad\\L2\\Vd\\L6 

<\\Vu\\l2{\\Vu\\l2 + \\VMl^)\\M\l^ (by (mH) with r = 6) (^.24) 
<e\\VMh + C{e)\\Vu\\l4Ad\\l, + \\Vu\\l, + \\Ad\\l,. 

To estimate the last term in the right hand side of (|2.18p . denote M{d) = Vd Vd - ^\Vd\H3. 
Then we have (Ad, Vd) = div(M(d)) and 

- Uf {Ad, Vd) dx = - Uf {{M{d), v)) da + M{d) : Vut dx 
J Jan J 

M{d) : Vut dx (since ut ■ v = -q- = on dVl) 
=^ j M{d) :Vudx- J {M{d))t : Vudx (2.25) 
<^ J M{d) :Vudx + C J |Vdt||Vd||Vn|dx 
<^ J M{d) :Vudx + Ci€)\\Vd\\l^\\Vu\\l2+€\\Vdt\\l2. 
Putting ([2:20]) - ([2:25] ) into (l2T8]l . we obtain 

\li J ^ + + + J P\u\^dx 

<^ J (M(d) : Vn + P(p)divtx) dx + e{\\Vdt\\l2 + \\V\\\h) (2.26) 

+ C {mu)\\L^ + llVnlli^ + llVdIlioo + 1) \\Vu\\l, 
+ C\\Vu\\l, \\Vp\\l, + C(l + ||V^||i.)||Ad||i. + C(6). 
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Now we want to estimate ||V/3||^2- Differentiating the equation (jl.ip with respect to x, multi- 
plying the resulting equation by 2V/3 and integrating over il, we obtain 

• Vp) dx 

= - j \Vp\^dbrudx - 2 j V{u) : V p (^V pdx - 2 J pVdivu-Vpdx '^^'^''^ 

<(\\V{u)U^ + l)\\Vp\\l,+e\\VMh- 

Next we want to estimate ||Vdf||^2- To do this, we multiply (j2.15p by Vdt and integrate over 
and use ^ = on d^l to obtain 

^ j \Ad\'^dx + j \Vdt\'^dx = j {V{\Vd\'^d) -V{u-Vd))Vdtdx 

<C{e) J (|Vu|2|Vd|2 + \u\^\V^d\^ + \Vdf + \Vd\^\V^d\^) dx + e||Vdt|||2 (2.28) 

<e||Vd4||i2 + C + C7||V(i||i..(||V2d||2, + ||Vn||i2)+C / \u\^\V^d\^dx, 



where we have used ()2.14p (with r = 6) in the last step. For the last term in the right hand side of 
(|2.28p . by using Nirenberg's interpolation inequality and (|2.14p we have 

\u\^\VM^dx <||n||i6||v2d||i3 < \\Vu\\UVdh4'^^dh2 + \\Vn\\l, 

<l|Vn||i2||V3d||^. + llV^lli^ < C{e)\\Vn\\l,+e\\Vmi,+C. 

Applying the standard //^-estimate to the Neumann boundary value problem of the equation ()2.15p . 
and using (|2.14p . we have 

llV^dlli^ <||VAd||i2 + llVdll^. < WVdtWl, + ||V(n • Vd)\\l, + \\V{\Vd\'d)\\l, + ||Vd||^, 
<\\Vdt\\l. + \\Vd\\l^{\\Vu\\l, + \\Vmi,) 

+ llVdll^e + ||Vd|||i + j \u\^\V^d\^dx (2-30) 
<\\^dt\\h + llVdIlioo (||V.x||i2 + \\V'd\\l,) + I \u\'\V'df dx + llV^dlli^ + 1. 
Substituting (|2.29p into (j2.30p and choosing e sufficiently small, we have 

llV^dlli. < 1 + ||Vd,||i2 + llVdIlioo {\\Vu\\l, + \\V'd\\l.) + \\Vu\\l, + llV^dlli^. (2.31) 
Substituting (fOTT) into ([2:29]l . we obtain 

J \u\'\V^d\' <C + e\\Vdt\\h + C|l Vdllioo (llV^zlli^ + II V^dlli^) + C||Vtx||i2. (2.32) 
Putting (|2.32p into (|2.28p and choosing e sufficiently small, we obtain 

1 1 \Ad\'dx + J \Vdt\'dx<l + \\Vu\\l, + \\Vd\\l^{\\V'd\\l, + \\Vu\\l,) + \^^^^^^ 

<1 + \\Vu\\l + \\Vd\\l^ {\\Vd\\l, + \\Ad\\l, + \\Vu\\l) (2.33) 
+ l|Ad||i2. 
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Putting (|2:26D . (12:2711 and (12331) together, we obtain 

(/i|V X'u|2 + (2^ + A)|divn|2 + |Vp|2 + |A(i|2) dx + J {2p\u\'^ + \Vdt\'^) dx 

<2- / (M(d) : Vn + P(p)divn) + e||Vdt||i2 + e||V\||i2 

+ llV^xlli^ + mu)U^)\\Vp\\l, +C{1 + \\V{u)\\l^ + llVnlli^ + ||Vd||i..) \\Vu\\l, 

+c{i + \\vu\\l, + ||vd||i..)||Ad||i2 + c\\vd\\l^ + a 

By VF'^'^-estimate of the Lame equation under the Navier-shp boundary condition ()1.10p (see [15 
Lemma 2.3 and also the proof of Lemma 2.2), we obtain, by using the equation (|1.2p and ([27 



<\\C^\\h + llVnlli^ < llVnlli^ + WpiiWh + l|V(P(p))||i2 + \\Ad. Vd\\l, 
<||Vn||i2 + \\pn\\h + llVplli^ + \\Vd\\U\Ad\\l,. 

Choosing sufficiently small e > 0, we have 

j {p\V X up + {2p + A)|divn|2 + |Vp|2 + jAdp) dx + j (2/)|n|2 + \Vdt\'^) dx 
<2— / (M(d) : Vn - P{p)divu) dx 

+C [1 + mu)\\L^ + llVnlli^ + llVdIlio.] [\\Vu\\l, + llVplli^ + IIAdlli^] + C{1 + ||Vd||ioo). 
Integrating from to t, < t < X,, and applying ()2.1ip and (12. 7p . we obtain 

(|Vu|2 + |Vp|2 + |Ad|2) r /" (p|n|2 + |Vdt|2) dxds 

(|M(d)||Vn| + |P(/5)||divn|)(t)dx+ / (|M(do)||Vno| + |P(/9o)||divno|) dx 



< 



+ y (iV-uol^ + iVpol^ + lAdoP) ^^ + '^_^ (1 + ||Vd||i^)ds 

+ /*[! + ||P(n)||L^ + llVnlli^ + llVdIlioo] [\\Vu\\l, + \\Vp\\l, + \\Ad\\l,] ds 
Jo 

<C+'^\\Vu\\l,+C l{\Vd\' + \p\)dx 

+C f [1 + \\V{u)\\l^ + llVnlli^ + llVdlli^] [llVnlli^ + llVplli^ + IIAdlli^] ds. 



(2.35) 



Since the coefficient function 

[1 + P(n)||L- + llVnlli^ + llVdlli^] G Li([0,r,]), 
the Gronwall's inequality. Lemma 2.3 and the conservation of mass imply that for any < t < T^,, 

/"(|Vn|2 + |V/9|2 + |v2d|2)(t)(ix+ / I {p\u\'^ + \Vdt\'^)dxds<C. 
J Jo Jn 

The proof is complete. □ 
As an immediate consequence of the proof of Lemma 12.41 we have 
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Corollary 2.5 Under the same assumptions as in Lemma \2.4\ we have 

sup \\dt\\l2 + [ * {\\pt\\l2 + \\Vu\\l^ + \\Vd\\j^2) dt < C. (2.36) 

0<t<T, Jo 

Proof. It follows from Lemma [2^ that Vu, V/>, Ad G LfLl{Q, x [0, T*]). By Sobolev's inequality, 
we then have u G LfLl{Vl x [0, T*]). On the other hand, Lemma O implies Vd G LfLl{n x [0, T^]) 
for r = 3, 4. Therefore, by (jl.3p . we have 



< (l^llVdl + |Ad| + iVdl^) G L~L^(J] X [0,r,]). 

It is easy to see that L^L^-estimate of V^u and V^d follows from (I2.31j) , (I2.34p . (2.1), and Lemma 
2.4i To see LfL^-estimate of pt, note that (jl.ip and Lemma |2 . 1 1 imply 



\pt\ < |Vp||n| +p|divn| < \u\\Vp\ + \Vu\. 
By the Sobolev's embedding, we have u G LlLf{Q. x [0,r*]). Hence 

lll''^l|Vp|||L2(nx[0,T,]) < ll^llL2Lg°(nx[0,T,])II^PllL?°L2(nx[0,T,]) < C. 

This clearly implies ||L2(r2x[o,T,]) ^ C. The proof is complete. □ 

Step 5. Estimates of {^ut,V^u,Vdt,V^d) in LfLl{Q. x [0,r*]). More precisely we have 

Lemma 2.6 Let < T* < +oo he the maximum time for a strong solution {p,u,d) to ^.l\) - fT7^) . 
[n\) together with fig) or or [TTTOl) . // [TTl]) and hold, then 



sup [ (p|utp + |V\p + iVdtp + iV^dp) dx + [ ' [ (iVujp + |djtp)dx(it < C. (2.37) 



Proof. For simplicity, we only consider the Navier-slip boundary condition p.lOp . Differentiating 
the equation ()1.2p with respect to t, we get 

putt + Ptut + pu ■ Vut + put ■ Vu + ptu ■ Vu + V{{P{p))t) 

(2.oo) 

= Cut-V ■ {Vdt (S)Vd + Vd(S) Vdt - Vd ■ Vdth) . 
Since ut ■ v = ^ and curlu^ x i/ = on as in the proof of Lemma 12.41 we can verify 

{Cut,ut) dx = y {p\V X ut\^ + {2p + A)|divutp) dx. 

Thus, multiplying (j2.38p by ut and integrating the resulting equation over O and using (jl.ip . we 
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obtain, by using Sobolev's inequality, Holder's inequality, and (j2.22p . 

~ J p\ut\'^dx + J (^|V X ntp + (2^ + A)|divnt|2) 

{p\u\\Vut\\ut\+ p\u\\V{u-Vu-ut)\ + \{P{p))t\\dwut\) dx 
+ / p\ut\'^\Vu\dx+ / |Vdt||Vd||Vnt|(ix 



< 



<\\Vut\\L2\\^Ut\\L2\\u\\L^ + \\Vu\\LA\ut\\LA\^nt\\L^ + \\V Ut\\ l2\\V dt\\ d\\ Lo 

+ ||pt|U2||divnt||i2 + / {p\u\\Vu\^\ut\+ p\u\'^\V^u\\ut\+ p\u\'^\Vu\\Vut\)dx 

2 

L6 



<\\^Ut\\LA\y/pnt\\LA\^u\\Hi + \\^u\\LA\^Ut\\L2\\u\\i, + \\ut\\ lA\^ ' u\\ l2\\u 

+ WutUAMLAWAh + WptUAW^MU^ + \\Vut\\LA\^dt\\LAWd\\L^ 

<l|Vnt||i2 (IIVp^dlL^llVnllHi + WVuWur + WptWv^ + \\V dtUA^^ d^^) 

1 



< 
~2 



p\yut\''dx + \\^put\\lA\Vu\\l, + \\vu\\l, + WptWl, + \\vdt\\U\yd\\l^. 



This gives 

^ j p\ut\'^dx + j (/z|V X + (2/i + A)|divnt| 
<||Vn||2,, / p|^,|2d^ + ||v2n||2, + ||p,||2^ + ||Vd||ioo||Vdt||i. + l. 



dx 



(2.39) 



Differentiating the equation (11.3P with respect to t, multiplying (it^ and integrating over J7, we 
dv 



obtain, by using ^ = on 517, Sobolev and Holder inequalities. Lemma 12.31 ^-iid Lemma [27 



\jt j 1^^*!'^^ + / l^«l''^^= J{dt{\Vd\^d-u-Vd) ,du)dx 

SldtthAWthAl'^dh, + WdtthAHlLAl^dth^ 
+ {i + \\vd\\U)\\vdt\\lA, 



<\\\du\\h + C[\\Vut\\l, + \\Vdth2\\V%h2 



which implies 

^ y \vdt\^dx+ J \du\^dx 



(2.40) 



<^^[||Vnt||i2 + \\V'dth2\\Vdth2 + (1 + ||Vd||ioo)||VdtrL2]. 

Now we need to estimate ||V^(it||i2. Li fact, by applying the standard //^-estimate on the 
equation (|1.3p and Lemma [2731 we have 

\\V^dt\\L2 <\Ndt\\L2 + Wduh^ + Wdtiu ■ Vd)\\L2 + \\dti\Vd\^d)h2 

<\\^dt\\L2 + Wdtth^ + WuthAl^dUs + \\u\\L4Vdt\\L^ 
+ \\dt\\L4'^d\\l, + \\VdthANd\\Le 

<\\du\\L2 + \\Vut\\L2 + WVdtWlAlVdtWl + \\vdth2 

<l\\^^dt\\L2 + C [\\du\\L2 + \\^Ut\\L2 + WVdtWLA . 
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Thus 

72 



W"^ dt\\L2 <\\dtt\\L^ + ||Vnt||i2 + \\Vdt\\L2. (2.41) 
Substituting ()2.4ip into (j2.40p . and using Cauchy inequality, we obtain 

I J \Vdtfdx + I |dtipdx<i||d«||2,+C7[||Vni||2, + (l + ||Vd||ioo)||Vdi||i.]. 



Thus 

d 
di 



J \Vdt\^dx + J \dtt\^dx <C [\\Vut\\l2 + (1 + \\yd\\loo)\\Vdt\\l2\ . (2.42) 

Multiplying (j2.42p by ^ and adding the resulting inequahty into (j2.39p , applying Lemma 12.31 and 
Lemma [231 and then employing Gronwall's inequality and applying (|2.1ip (with u replaced by Ut), 
we obtain ^ 

sup I {p\ut\^ + \Vdt\^)dx+ i {\Vut\^ + \du\'^)dxdt <C. 
o<t<T, J Jo Jn 

To estimate V'^d in L^°L'^{Q x [0,T^]), first observe that by Nirenberg's interpolation inequality, 



we have 



\\vdU^<\\vdh2 + \\vd\\Uv^\l,. 



L2\\ ^ "IIL2- 

Putting this inequality into (j2.3ip and using L^L^-bounds of Vdt, Vu, V'^d, we obtain that for any 

< t < r*, 

\\v^d\\l,<c + c\\v^d\\l<^\\v'd\\l,+c, 

which clearly yields that 

sup (iivdiUoo + ||v3d||i2) < a 

0<t<T* 

To see V'^u G Lf°L'^{Q, x [0,r,,]), observe that the ff^-estimate on the equation (|1.2p under 
(fTTO]) . ([27711 . and Lemma [23] imply that for any < t < T*, 

I|V^u||l2 ^l|Vn||i2 + \\Cu\\l2 

<1 + \\,/put\\L2 + ||n • Vu\\l2 + \\V{P{p))\\l2 + \\V^d\\L2\\Vd\\L^ 

1 1 

<1 + \\^/pUt\\L2 + ||^t||L6||Vu||i3 < 1 + \\^/pUt\\L2 + || Vn|| ^2 || V^ujl^a 



<^||V\||i2 +C. 



1 

^2' 

In particular, we have 

sup ||V^m||2,2 < C. 

0<t<Tt. 

The proof is now complete. □ 

Step 6. Estimate of Vp in Lf°L%{Q x [0,T,,]) for some 3 < g < 6. With the estimates already 
established by the previous Lemmas, we then have the following Lemma. 

Lemma 2.7 Let < < +oo be the maximum time for a strong solution {p,u,d) to il.l\) - [T73\) . 
together with fig) or or (TTTOl) . // [Til]) and fOP hold, then 

sup (maxWptUr + \\p\\wi>-i] + f {\\u\\D^,, + \\^^dt\\l2 + \\V^d\\l2)dt<C, (2.43) 
for any 3 < q < 6. 
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Proof. For 3 < g < 6, by the same calculations as in [15] Lemma 2.5, we have 

(|Vy9|«)t + div (iVpl^u) + {q- l)|Vp|'?divu + q\V pl^-^iV pYV{u){V p) + qplVpl^-^Vp ■ Vdiv?/ = 0, 

which yields that for ([L8]) or ([L9D 



^IIVpIIl, < C {\\V{u)\\loo + 1) IIVpIIl, + C||Vdiv^/||M, (2.44) 
at 

and that for (fTTOl) 

|||Vp|U, < C {\\V{u)\\l^ + 1) IIVpIIl, + C||VG|U,, (2.45) 

where G = (2/i + A)divii — P{p). 

For boundary conditions (jl.Sp or ()1.9p . by using the L^-estimate for the elliptic equation and 
(|2.7|) we have 

llV^nlli, <\\put\\L^ + ll/o-u • Vu||l, + ||V(P(p))||m + ||Ad • Vd^i + 1 

<II^/P^dl^lkdiy + hllL-l|Vu||i, + IIVpIIl. + ||Vd||ioo||Ad||i, + 1 (2-46) 

<\\Vut\\L^ + \\Vp\\L', + l. 

Substituting (I2.46P into (I2.44[) . and using Gronwall's inequality, we obtain the bound sup ||/o||iyi''? 

0<t<T* 

for the first two boundary conditions (jl.Sp and ()1.9p . 
For boundary condition (jl.lOp . we rewrite ()1.2p as 

VG = pV X curln + put + pu-Vu + Ad - Vd, (2.47) 

which yields that G satisfies 

I AG = div {put + pu-Vu + Vd- Ad) , in Q, 
1 VG • ly = —p{u ■ V)i^ • n, on dQ, 

where we have used that (V x curlu) • vIqq = (curln x uIqq = implies (V x curlu) • vIqq = 0, 
see [15] page 33 or [5l [6]), Vd ■ i^ldn = and u ■ iy\dn = 0. 

Using the L^-estimate for Neumann problem to the elliptic equation (I2.48p . we have 

\\'^G\\l'i Slputh'i + \\pu-Vu\\Lq + \\Vd- Adlli'i + \\p\u\'^\\c(n) 

<\\^uth^ + l. ^'^^ 

Puting (|2.49p into (|2.45p . we obtain the bound sup ||/o||vf1''? by Gronwall's inequality. 

0<i<T, 

For r = 2 or (jl.ip implies 

\\Pt\\L^ ^ ||tt||L°°||Vp||L'- + ||/0||L°°||divii||Lr 

< II Vnll^i llVplUr + WpWl^. \\Vu\\hi < C. 

It follows from ^^AT\i that 

\\'^^dt\\l2<\\dtt\\l2 + \\Vut\\l2 + l. 
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This, with the help of (|2.42p . imphes, after integrating over [0, T,,,], 



[ \\V^dt\\l2dt<C. 
Jo 



(2.50) 



Applying the standard L^-estimate to ()1.3p . we have 

llV^dlli. < \\V%\\l, + \\V\u-Vd)\\l, + \\V\\Vd\'d)\\l, + l 

< WV^dtWl^ + \\u\\loo\\V^d\\l2 + \\Vd\\la.\\V^u\\l2 + \\Vu\\l4V^d\\l3 + 1 

< \\V^dt\\l2+l. 

Integrating this inequality over [0,T^,], and using (|2.50p . we get 

\VU\\l2dt < C. 

'0 

By the bound on \\V p\\li in ([233]), and (fTHT]) . we easily see that 

(•T, 



J 

Jo 



[ ' \\V^u\\l,dt < C. 
Jo 



holds for (jl.Sp or (jl.9p . For the boundary condition (jl.lOp . since u • = on dfl, it follows from 
Bourguignon-Brezis [2\ (see also [15j Lemma 2.3) and (|2.7p that 

< ||V(divn)||L9 + ||V(curln)||L, + ||Vu||l5 

< IIVGIIl, + ||V/5||m + \\Vu\\h^ + ||V(curltx)||i, 

< l + ||VG||M + ||V(curl?x)||L,. 

Since (V x uY = on dQ, it follows from [33j that 

||V(curln)||j;,<j < ||div (curl n) 11^9 + II V X curl n||ig < II V X curl u||i<7, 
where we have used the fact that div (curl u) = 0. On the other hand, since 

X curl u = VG — puf — pu ■ Vu — Ad ■ Vd, 

([Tigj) implies 

||V X curlnlli, < 1 + ||Vnt||i2. 
Putting these estimates together, we have 

I|V^u||l. <l + \\Vut\\L^ + ||VG||l,, 

which clearly implies 



|v\||i, < c. 







The proof is now complete. □ 

Step 7. Completion of proof of Theorem 11.31 

With the above established estimates, we obtain (2.2) and (2.3). This implies that is not the 

maximum time of existence of strong solutions, which contradicts the definition of T^,. Therefore, 

(j2.ip is false. The proof of Theorem 11.31 is now complete. □ 
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